General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It wa 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or picti 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of s 
of the material. However, it is the best reproduction available from the originc 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




^ II8TITIT1 IE PEigUISIS ESHCIMS 


6. Palavras chaves - selecionadas pe1o(s) autor(es 
DIFUSHO; IMAOEm LAWSAT; POLVIffXO 

) 

7. C.D.U.: SSI. 511. 6:551.510.42 
r Z77-: 



1, publica^ao n9 

XNPE-’BSSS’-PRE/OdB 


2. Versao 


3. Data 

Fev,,19B2 


4. Origern 

DME 


Programa 

RECAT 


5. Distrlbulgao 

□ Interna 03 Externa 

□ Restrita 


ESTUDO DA DIFVSAO ATMOSFFRICA 
VSARDO IMAGERY DO SATELITE LAWSAT 


9* Autoria Joae Apaveoido Torsani 
Yeliasety Viawanadham 


'll. Oltima pagina: 29 


12. Revisada por 



^avloa de ''Almeida 


13. Autorizada por 


Assinatura responsavel 

14. Resumo/Notas 

0 modelo gauaaiano e aomumente utilizado para aaloular a 
oonaentragdo de jpoluentea. Vara iaao ado neoeaadrioa oa valorea^ doa para 
metroa de difuaao (deavio padrdo da diatribuigdo de oonaentragdo) para 
diferentea oondigoea, Neste trahalho, valorea do dea'Oio padrdo da diatri^ 
buigdo da oonaentragdo , ao longo do eixo Y favam medidoa a partir 
de imagena obtidaa pelo aatelite LANDS AT e ^ oomparados oom os valorea 
previatoa peloa modeloa de Paaquill^Giffordf Julioh e Hogatrdm. 





Nelaorr de Jeaua Parada 


Dirotor 


SHOR U.wsmi58!fs 


©M5®. gSEiTftS 




15. Observagoes Submetido d Reviata Cimoia e Cultura 










ESTUDO DA DIFUSAO ATMOSFERICA 
USANDO IMA6ENS DO SATELITE LANDSAT 


J.A. Torsani^^) e Y. Visv^anadham^*) 

Conselho Nacional de Des?Dnvolvimento CientTfico e Tecnologico - CNPq 
Institute de Pesquisas Espaciais - INPF 
12.200 « Sao Jose dos Cavn;^os - SP - Brasil 


Resu mo 

0 model 0 gauss iano e comumente utillzado para calcular a 
coneentragao de poluentes. Para isso sao necessarios os valores dos p£ 
ranietros de difusao (desvio padrao da distribuigao de concentragao) p^ 
ra diferentes condigoes. Neste trabalho, valores do desvio padrao da 
distribuigao da concentragao, ao longo do eixo Y (oy)» foram medidos a 
parti r de imagens obtidas pelo satelite LANDSAT, e comparados com os 
lores previstos pelos modelos de Pasqui 11 -Gifford, Julich e Hogstrbmi 


(1) Assistente de Pesquisa, INPE, S.J. Campos, S.P. 

( 2 ) Pesquisador Pleno, INPE, S.J. Campos, S.P. 



1 . Introdugao 


Os poluentes sap introduztdos no meio amblente atraves 
de diversas atividades, causando efeitos de larga escala e amplo espa 

t *** 

lhainento, Incluindo os efeitos sobre o cllma da Terra; os efeitos de pe 
quena escala estao assocfados a poluigiTo urbana, cujo transporte e dj[ 
lulgao sao devidos aos niovinientos do ar. 

Em anos mais recentes, os metodos para estimativa de dis 
persao atniosferfca forani modiftcados, com base em resultados obtidos 
por medidas experimentais, Pasquill (1961) propos urn metodo para a est^ 
mativa de dispersao. Este metodo, com algumas conversoes feitas por 
Gifford (1961), e atualmente amplamente utillzado, Em vista dos proble 
mas que envoi vein a redugao dos dados meteorologicos para o calculo dos 
parametros de difusao, tem-se procurado determinar tals parametros di^ 
retamente das plumas, utilizando-se tecnicas de sensoriamento remote. 

0 uso do sensoriamento renioto para monitorar e identif[ 
car fontes poluidoras tern sido explorado nos ultimos anos. Randerson 
(1968) publicou futografias, obtidas de satelites tri pul ados, que mo^ 
tram plumas visTveis perto da costa do Texas, 

Com 0 uso de urn analisador automatico (GE Image TOO) e 
de imagens obtidas do satellte LANDSAT, foi possTvel obter perfis de 
plumas com extensao de* aproximadamente 130 km. Desta forma, imagens do 
LANDSAT podem ser usadas para detectar e monitorar plumas de grande ex 
tensao, 

0 objetivo deste trabalho e obter medidas do desvio p^ 
drao da cHstribuigao de concentragao com relagao ao eixo y (Oy), a par 
tir de imagens obtidas pelo satelite LANDSAT. Tendo-se os valores de 
Oy*, obtidos das fotografias, comparar-se-‘ao tais valores com os previ^ 
tos pelos modelos de Pasquill -Gifford, JUlich e Hdgstrdm. 
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2. TeoHa Basica 

A apltca 9 ao pr^tlca da equa^w da difusao requer o conhe 
cimento dos desvios padi'oes Oy e 0^1 Esses parametros sao dados em fuji 
9 S 0 da distancia da fonte, das condigoes de establlidade da atmosfera, 
da aVtura da emissao da fonte e da nugosidade-superficlal . Alguns dos 

1 

mals importantes sistemas Incorporados na apUcagao pnatica do compute 
da difusao sao: 

2.1 - Sistema de Pasqui 11 -Giffo rd 

Pasquill 0^61 J classificou seis tipos de difusao, depeji 
dendo da establlidade, as quais podem ser determinadas pela velocidade 
do vento - medida a uma altura de 10 m, e pela insolagao - definida para 
perTodos do dia e grau de cobertura das nuvens, A expressao proposta 
para 0 calculo de Oy, em fungao da distancia, e dada por: 

ffy (X) = (ai Inx + azi x Cl) 

onde os coeficientes sao especificados na label a 1 . 

2.2 - Sistema JUlich 


As formulagoes apresentadas neste sistema foram elab^ 
radas para nTveis de emissao de 50 e 100 m, Neste caso, 0 nivel de emi£ 
sao e de 76 m; por isto houve uma interpolagao linear, cc>nforme propo£ 
ta feita por Vogt (1977). De inTcIo a classificagao foi organizada de 

tal forma que uma comparagao com as categorias de establlidade de 

Pasquill -Gifford fosse possTvel . 0 parametro a e dado pela fungao 

3y (X) = Py X ‘>y 

onde Py e cjy sao coeficientes e estao listados na Tabela 1 . 
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ORlGd^AL ITi 
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Tabela 1. Coeficientes de difusao para os dlferentes sistemas e catego 
rias de estatiilidade* 


SISTEMA DE 

COEFICIENTES 


CATEGORIAS DE ESTADIUnADE 1 

DIFUSAO 

DE DIFUSAO 

A 

8 

1 C 

0 

£ 

F 

Pasquill- 


-0,0234 

-0,0147 


-0,0059 

-o.oqs'i 

-0,0029 

Gifford 


0,3500 

0,2430 


0,1030 

0,0330 _ 

0,0540 




0,4094 

0,3534 


1 ,69T0 

5,3320 

jo1irf> - 

y 

■■■■ 



■HB 



V«lor«s In 








tcrpolados' 








para 76m 

‘’y 

B 

3,8991 

0,8733 


0,6211 

0,5778 ■ 


,9321 


3,8991 


0,8733 


0,621 


0,5778 ■ 
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2.3 - Sistema Ndgstram 

Neste sistema o calculo de oy efetivo e fetto usando-se 
a seguinte expressao (Hdgstrbm, 19641: 

Oy »Oy^ + (3) 

pnde: 

a^d “ variancia devida 3 difusao tur&ulenta 

o^s s variancia devida ao cisalhamento vertical da diregao do 
vento (este fator e signiftcativo para longas distancias) , 

cTyd e Oyg sao obtidos a parti r das seguintes relagoes: 

Cyd (X) = 122 S /2 {exp t-Coxl + Cox - 1 D" t^) 

onde : 

Co = 10"^ 

'o,8 + para condigoes estaveis e.neutras 

g ^ J 1 f 3,4 . 10“^ Rg 

1,4 para condigoes instaveis 

onde R^ eoparametro de estabilidade dado por: 
s * 

R . ( J!I + r ) . 10* C5) 

5 3Z 

sendo que: 

•51 = Gradiente vertical de temperature. 
dZ 

r' ■ Taxa de variagao vertical da temperatura adiabatica seca. 

« Vento medio na altura de 500 metros. 


0 termo devido ao cisalhamento e dado por: 
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(x) dx - 



( 6 ) 


onde: 

c = 1 ,Q para a- ^ — 

^ 2 

9^7 9 h 

s -i- + -£- para — a, <: 2h 

3 h 3 ’ 2 ^ 

= 2 para ^ 2h 

\p = cisalhamento vertical da diregao do vento (rad. m’^). 

0 (x) a desvio padrao da distribuigao de concentragao ao iongo do 
eixo vertical. 


Para o caso neutro, o desvio padrao vertical e expresso pel a formula: 


(x) a 30,4 [2 (exp (-ao x) -f ao x 



( 7 ) 


onde: 


ao = [131,02 log 

V 

Neste trabalho, h. = 76m e 7^ = lO'^m (CounIham, 1975). 
ra 0 caso estavel , pode-se escrever o desvio padrao vertical como: 



zn 


(1 + a R ) 


( 8 ) 


Na Equagao8,R eo parametro de estabilidade, ja definido.A, ,quaD.tidade..a 
e dada por: 


a a 6.10-®.h^«P'^ .Ca“^ 


onde: 

Ca a coeficiente de arrastamento geostrofico, referente a condj^ 
goes adlabaticas, dada pela formula: 



• 6 " 


log (Cg Ro^j) -2,24 

On X *9»9 

® Zo.f 

onde Ro^ e o Numero de Rosshy superficial j Vg»o e a velocidade do vento 
geostrofico na super, fTcie; f e o parametrodeCo*r.i o1 is (« 2 n seii(() s*^); 
SI e a frequencia angular de rotagao da Terra; e (|) e a latitude. 

A forn^ulagao para o caso instavel nao sera apresentada 
neste trahalho, por nao se ter encontrado nenhum caso para esta condj[ 
gao de estabil idade. 

3. Determinagao da Categoria de Estabil idade Atmosferica 

0 metodo mais comuinente utilizado e aquele proposto por 
Pasquill^uifford (1961), onde a categoria de estabil idade e determinada 
utilizando-se: 1) Velocidade do vento no nTvel de 10m e 2) Estimativa 
quantitativa da insol agao durante o dia, e cobertura de nuvens durante 
a noite (noite refere-se ao perTodo de uma hora antes do por do sol e 
uma hora depois do nascer do sol. 

As condigoes de estabil idade sao divididas em seis cate 
gorias, conforme Turner (1970): 

A “ Extremamente instavel. 

•B - Moderadamente instavel. 

C ~ Levemente instavel. 

D “ Neutra. 

E • Levemente estavel. 

F - Moderadamente estavel . 

4. Local de Estudo Dados DisponTveis 


0 local de estudo foi Cabo Frio, local izado na latitude 
de 230 00‘ S e longitude 42o 00’ W. Este local e referente a orbita 



m ponto 28 do satelite LANDSAT. 
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Nas Imagens vtsuallza-se uma plutna provenlente de uma 
fonte puntual (Companhia Nactonal de AlcaHsJ, situada no Arralal do 
Cafao. Os parametros de emissaO} conforme Informagao prestada pel a Com 
panhia, sao: 


/ ^ 

Altura da chamlne 

Vazao dos gases ^ . 

Velocldade de saTda dos gases . 

Composlgao dos gases em volume: 

Vapor d^agua 

N2 

Oa 

CO2 


76 m 

20Q m* S-* 
7 m s^^ 

47^ 

m 

m 

2% 


CO, NO^, SO2, Mg, MgSOn, NaCl Tragos 


Devido a necessidade de oUter perfis verticals da velo 
cidade do vento e da temperatura, e tendo-se em vista que no local de 
estudo conta-se apenas com os dados de superfTcie, utilizam-se radios 
sondagens do aeroporto do Galeao, local izado na latitude de 229 5Ct' S e 
longitude 430 14' W, Rio de Janeiro ^ RJ, ohtidas no horario das 11:00 
e das 23:00 GMT, Cabo Frio encontra-se aproximadamente a 100 km Teste 
do Rio de Janeiro (Figura 1)^. 

Os dados necessarios § realizagao dos calculos estao 
descritos a seguir: 


0 primeiro con junto de dados utilizados esta relacionado com ob 
servagoes* de superfTcie, realizados na estagao meteorologica da Comp^ 
nhia Nacional de Alcalis, assim como os dados obtidos de radiossond^ 
gens realizadas no aeroporto do Galeao, no horario das 11:00 GMT, e os 
dados de cobertura de nuvens (porcentagem) e elevagiio solar, fornecidos 
pelo satelite LANDSAT, 0 conjunto de dados citados acima encontra-se di^ 
posto na Tabela 2. 




OCEANO ATLANTICO 



T 
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TabeVa 2. Parametros utlllzados para determlnar as condigoes de esta 
billdade para os tree casos escolhidos, 



DATA 

TEHPERATURA (oc) 

T.V.ViT 

(OC/lOOm) 

VERTO SUPERFICIAt 



COBERTURA DE 
NUVEHS {%) 

DO AR 
(T.) 

OA AGUA SUPERFICIAl 
(Th) 

OIREQAO 

VEIOCIOAOE 
(ms* ‘ ) 

01 

18/4/75 

25,0 

22,5 

-0,72 

N 

1,0 

5,8 

360 

20 

02 

23/6/7C 

23,0 

22,0 

0,56 

NNE 

3,3 

4,8 

220 


03 

22/5/77 

21,0 

O 

CVI 

0,47 

W 

1,1 

4,7 

300 

40 


Uf « velocidado do yonto tornado no nTvcl dc 500m . 
T.V.V.Ti » Tdxa de varia^ao vertical da tcmporatura. 
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0 segundo grupo de dados compreende valoros de tempera 
tura da agua superficial do mar, velocidade e diregao do vento para pe 
rTodos de uma semanai na regiao de Cabo Frio, os quais saoapreser.tados 
na label a 3, 


5* Imagens LANDSAT 

Os sate! lies da seric LANDSAT anteriomente denomlnada 
ERTS, sao platafomias de sensoriamento remote de dados colocadas em 
orlxita terrestre de grande fncTinagao com o equador, quase-circular, a 
uma altitude de aproximadamente 918 km, A trajetoria orbital e feita 
de maneira sTncrona com o sol, sendo o perTodo orbital de 103 min 16s, 
Este efeito faz com que, apos 251 revolugoes (isto e, a cada 18 dias), 
0 satelite volte a imagear o mesmo local. 


0 imageador multiespectrai MSS-LANDSAT, daqui para a 
frente denominado somente MSS, e urn dos subsistemas do satelite e con 
siste num espelho oscilatorio e num sistema o'tico que reflete a radi% 
^'ti{ r:nda da cena, num conjunto de 24 detetores divididos em 4 bandas 
o« canals, cada urn com 6 detetores, Cada banda e sensTvel a uma deter 
minada faixa do espectro, a saber: 


canal 4 - 0,5 a 0,6 ym 1 
canal 5 - 0,6 a Q,7 pm J visTvel 


canal 6 
canal 7 


0,7 a 0,8 ym I -j^if^^a-vermelho proximo 
0,8 a 1 ,1 ym J 


5.1 - Sistema I -100 


0 s'istema interativo de analise de imagem multiespectrai 
(I-lOO) e urn sistema que realiza a classificagao semi-automatica de 
Imagens, extraindo informagoes tematicas de toda a cena ou de parte de 
la. A entrada da imagem no sistema se da por meio de fitas magneticas 
compatTveis com computadores (CCT), que fornecem: 

a) 0 nTvel de cirrza em cada canal do MSS, para cada "pixel", o]i 
de 0 sistema acessa os quatro canals ao mesmo tempo, e 
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Tabela 3. Parainetros meteorologicos da supei’fTcie. 
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Clw. |»V1 
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22,0 

23.0 

12.1 

- 
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VSlOCtOAOt 
00 HNIO 
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mun 

\mtn 
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U/04/f» 

c 

jy ’ 
w 

N 

Nl 

)y 

1 

N 

0,0 

f,o 

*t0 

1.0 

^.0 

1.0 

>•0 

?J.0 

24.0 
!4,0 

25.0 

22.0 

23.0 ■ 

21.0 

. 4,0 

1.0 

4.0 

1.0 
1.0 
0,0 
2,0 
1,0 

24.0 

21.0 
21,0 
2C.0 
21,0 

24.0 

n.o 

21.0 

11. 0 

21.1 

22,1 ’ 
21,0 
21,0 

23.0 

21.1 
23,0 

IW 

ly 

4 

NC 

ly 

IW 

Nt 

N 

1 

24.0 

22.0 
21.0 

13.0 

23.0 

21.0 
21.0 
12,0 

lii.i 

22.0 

21.0 

23.0 

21.1 
21,0 

23.0 

23.0 

1 

li/w/jil 

■ t, 

0,0 

21,0 

22,1 

It 

M 

21,0 

13,0 

Ifc 

4,3 

22,0 

23,0 

■ 


t 

0,0 

23.0 

23,0 


l.l 

21,0 

n,i 

41 

1,1 

22,0 

21,0 

D 


W 

l.l 

22,0 

22,1 

it 

2,1 

21,0 

21,0 

y ' 

2.0 

21,0 • 

23,0 

■ 
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M 
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22,1 

Nt 

2.1 
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22,1 

y 

M 

21,0 

21.1 

■ 
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l.l 
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mi 
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y 
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■ 
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Mht 
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hdt 

f.o 

24,0 

> 22.1 

y 

l.l 

21.0 

22,0 

1 
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NS 
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22,0 

■1 

l.l 
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22,0 


U/ov» 

n 

l.o 

2S,0 

23.1 

1 

7.1 

20,0 

J4.1 

■a 

mm 

22,0 

24,0 
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Hi 
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17,0 

21,1 

lY 

10,1 

21,0 

23,0 


Bw 

20,0 

23,1 

«> 
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■1 
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11,1 

s 
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IW 
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HI 
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1 

1.3 
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■SH 
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.1 
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11,0 

21,0 

It 

M 
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23,0 
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11,0 

23,0 
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M 
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20,9 

23,0 

Nt 

1.1 

21,0 

23.0 

HD 

3.3 

21,0 

23,0 

— 
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■1 

M 

21,0 

21,0 

1 

2,1 

2A,0 

21, 0‘ 


1,1 

21,0 

23,0 
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hi as coordenadas de cada “pixel** em relagao a matriz de 3240 x 
2340 *'pixels'* que compoem uma 1n)agen). 

A iniagenii depois de ter side acessada, pode ser projet^ 
da total ou parcialmente em uiu terminal de vTdeo colorido. Este’ term^ 
?^al funciona conio uma tela de televisao e possui: 512 x 512 pontos 
(“pixels" do I-^lOO), sendo que para cada ponto e transferido urn “pixel" 
da imagem. Dependendo do tamanho da area, que foi projetada no video, 
ocorrera uma ampliagao maior ou menor da imagem. 0 sistema realiza a 
classificagao seml-automatica das imagens, segundo urn metodo interati^ 
VO (anal ise/maquinal haseado em: 

al informagoes fornecidas pelo analista, as quais sao introduzi^ 
das por meio de areas de treinamento, ou seja, areas de cara£ 
teristicas multiespectrais conKecidas; 

hi programas de classificagao, atraves dos quais o analista pro 
cessa os dados introduzidos e procura alcangar os resultados 
desejados. 

0 resultado da classificagao aparece em tempo real no 
terminal de video colorido. Neste ponto, o analista interpreta o resuX, 
tado e, se quiser, pode modifica-lo atraves da utilizagao de urn outro 
programa, ate alcangar urn resultado satisfatorio. Ha uma interagao to 
tal. entre d analista e a .maquina, proporcionando uma real imentagao do 
sistema em urn curto espago de tempo. 

5.2 “ Processamento das Fotografias . 

As Imagens sao inicialmente digital izadas em fitas ma^ 
neticas e process/Jdas no I-IQO. Tendo-se em vista que cada imagem apre 
senta uma condigao espacTfica, escolheu~se por comparagao visual, o ca 
naloucanais com os melhores registros da pluma. Utilizou-se, tambem, 
urn dispositive de multipl icagao de tons de cinza a fim de roelhorar o 
contraste. 
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As imagens no l^lOiO sao fotografadas com f limes negat1_ 
VOS e/ou com diapositivos, dps quats se otitem ampllagoes em copias de 
papel. A partir das ampitagoes, estabelece-se o contorno aparente da 
pluma e fazem-se inedidas para a determinagao do parametro Oy em fungao 
da distancia da fonte, Assume-se que a dtstrihuigao de concentragao 
na pluma seja aproximadamente gauss i ana, ob.tendo-se entao, os valores 
de oy a par;tir da relagao Y(x) ® 2,15 Oy. onde Y(x) e a distancia hori^ 
zontal transversal medida do centro da pluma ((J. ) ate o contorno, onde 
se assume que a concentragao e aproximadamente 10% daquela do centro 
da pluma (X ^ 0,1 X Q. It conforme mostrado na Figura 2, 

6. Resultados 

Com a aplicagao do sensoriaroento remoto, utilizando-se 
imagens do satelite LANDSAT, real izaram-se medidas de Oy em tres dias 
diferentes* Nas Figuras 4, 6 e 8 estao plotados os valores de Oy* obtj^ 
dos atraves das fotografias, e os valores cal cul ados a parti r dos tres 
metodos propostos. Salienta-se que os metodos de Pasquill-6ifford(P-:G) 
e JUlich (J), utilizaram, para classificagao da categoria de estabili^ 
dade, os parametros metecrologicos apresentados por Turner (1970), (T^ 
bela 3.1), enquanto o metodo de Hdgstrdm utiliza o parametro definido 
na Fquagao 5, onde: 

1 < Rg < 30 para a categoria de estabilidade neutra; e 
Rg > 30 para a categoria de estabilidade estavel , 

Desta maneira, os tres casos serao discutidos do ponto 
de vista dos tres metodos. Calculou-se tambem a porcentagem de erro p^ 
ra as distancias de 1 , 5, 10, 15 e 20 km de distancia da fonte, e os 
resultados estao apresentados na label a 4. A seguir discutir-se-ao os 
resultados de cada caso. 

Caso n9 1 


Medidas de Oy toram feitas utilizando-se a fotografiaob 
tida do satelite LANDSAT do dia 18/4/75, Figura 3, e os valores enco]i 
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Figura 3 - Fotografia obtida da imagem do 
MSS do LANDSAT; Caso n9 1. 
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tr'am-se plotados na Figura 4. Para este dia a categorla de establlida 
de determinada pelo metodos de P-G fdt A/B. Entretanto, os valores me 
didos na fotografla apresentaram-se proximos daqueles previstas por 
P“G para a classe de estabilldade neutra (D), portanto, o metodo de 
P-G^ superestimou os valores de Oy (66^ para a distancia de 1 km, 163^ 
para 5 km e 136% para 10 kmi, conforme mostra a Tabela 4. 0 metodo de 
JUlicK, que utlHza a mesma categorla de estabilldade do metodo de P-^ 
•apresentou valores de Oy que superestimaram os valores medidos em 115%> 
265% e 246% para as distancias de 1, 5 e 10 km, respect ivamente. Ja o 
metodo de Hbgstrdm, que utiltza outro criterio para determlnar a cat£ 
gor.ia . de estabilldade, baseando-se na Taxa de Variagao Vertical d(i 
Temperatura (TVVTi e na Velocidade do vento no nTvel de 500 m C^quagao 
apresentou valores de Oy bastante concordantes com aqueles >nedj^ 
dos na fotografla, onde a porcentagem de erro obtida (p%, 28% e 3% p^ 
ra as distancias de 1, 5 e 10 km, respectivamente) e balxa. 

Embora se tenha feito uma comparagao dos valores de Oy 
com os metodos.descritos, ressalta-se que cada experimento foi realize 
do em condigoes diferentes. No experimento realizado por Pasquill 
(.1961), a emissao tomou lugar proximo ao nTvel do solo, com urn comprj^ 
mento de rugosidade em torno de 3 cm; e medidas de concentragao ao nj 
vel do solo foram feitas para distancias menores que 1 km da fonte. P^ 
ra .distancias ma lores que 1 km, foram feitas extrapolagoes baseadas em 
alguns dados experimental s. Quanto ao experimento realizado por Vogt e 
Geiss (1974) no Centro de Pesquisa Nuclear de JUlich, Republica Fede 
ral da Alemanha, para alturas de emissao de 50 e 100 metros, as medj^ 
das de concentragao foram feitas ate a distancia de 11 km da fonte, ori 
de parte do terreno era aravel e parte arborizada. Finalmente, o meto 
do de Hdgstrbm incorpora parametros importantes como a rugosidade siJ 
perficial , altura'de emissao da fonte, taxa de variagao vertical da 
temperatura, velocidade do vento acima da camada de fricgao, alem de 
considerar que o parametro Oy efetivo e obtido atraves de uma parcel a 
ajj, que e a variancia devida 3 difusao turfaulenta, e uma parcela 
Oy,^, que e devida ao cisalhamento vertical na diragao do vento e si£ 
nificativa para longas distancias da fonte, considerando-se ainda o 
efx^ito do arrastamento para condigoes estaveis. Com relagao aos baixos 
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Gifford, JUlich e Hdgstrdm. 
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Tabela 4. Porcentagem de erro (1) calcglada a partir dos 
valores inedidos na foto para os tres metodos. 


y 


CASO 

H9 

DATA 

VALOiiES DE oy 

PORCENTAGEM DE ERRO (S) 

.HCTO 

DO** 

OISTANCIA DA FONTE (km) 

mEto 

D0“ 

DSSTANCIA da FONTE (km) 

1 

5 

10 

15 

20 

1 

5 

10 

15 

20 



- ' ■—•■—'I 

FOTO. 

100 

260 

520 

- 

- 







01 

18/04/75 

HtiG. 

100 

332 

504 

- 

- 

hOo, 

0 

- 28 

3 

- 

• 



P-G 

166 

683 

1.230 

- 

- 

P-G 

- 66 

-163 

-136 

- 

* 



JUL. 

215 

950 

1.800 

- 

- 

jOl. 

-115 

-265 

-246 

- 

- 



FOTO 

72 

136 

210 

309 

454 







02 

23/0fi/76 

HOG, 

90 

300 

465 

611 

761 

hOg, 

- 25 

-118 

-113 

- 98 

- 68 



P-G 

67 

289 

537 

769 

991 

P-G 

7 

-ni 

-146 

-148 

-118 



jDl. 

110 

460 

827 

1.165 

1.487 

Jill.. 

- 64 

-236 

-279 

-277 

-227 



FOTO 

80 

401 


« 

. 







03 

22/05/77 

h5g, ^ 

90 

299 

456 

584 

706 

hOg. 

- 

- 

- 

- 

- 



P-G 

146 

614 

1.126 

1.600 

2.048 

P-G 

- 

- 


- 

- 



jUl. 

204 

876 

1.616 

2.327 

3.014 

jUl. 

- 

- 


- 

- 
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valores niedJdos na fotografda i deve-^se considerar: que se trata 
de difusao sobre unia superfTcie de agua. As caracterTstlcas de difusj^ 
vidade sobre a agua ainda sao pouco conhecldas, comparando*-se com aqua 
las sobre o contlnente. Neste trabalho os resultados foram obtidos pra^ 
ticamente sobre o oceanoj razao pel a qual os resultados previstos pelo 
nietodo de P-G sao bastante elevados* em comparagao com os valores rnedl^ 
dos atraves da fotografia, Resultados oNtIdos por Raynor et alii (1975) 
tambem mostrararn gue os valores de Oy obtidos de experimentos realiza 
dos sobre uma superfTcie de agua foram menores que os encontrados em e£ 
perimentos realizados sobre o continente, sendo que para condigoes in^ 
tave-'’- os valores estiveram abaixo dos prevtstos por P-Q para a catego 
ria e para condigoes estaveis os valores estiveram abaixo da catego 
ria F. Lyons (1975) verificou que, apos alcangarem a distancia de 30 a 
40 km sobre a superfTcie da agua, algumas plumas que originalmente tj[ 
nham uma caracterTstica instavsl (A), apresentavam caracterTsticas estS 
veis (F). 

Normal mente a intensidade de turbulencia sobre uma s^ 
perfTcie de agua e bem menor do que sobre o continente, devido: 

1) ao efeito de redugao da turbulencia mecanica, causada pela s£ 
perfTcie lisa da agua; 

2) as forgas de flutuagao resultantes da estabilizagao de massas 
de ar durante perTodos de advecgao de ar quente sobre agua 
fria; e 

3) ao grau de variagao da dispersao da pi uma ocorrido do continen 
te para a* agua ser dependente do comprimento alcangado pela 
pi uma sobre a regiao. 

Desta forma, metodos teoricos e experimentais sao nece^ 
sarios para explicar o fenomeno de dispersao na camada limite marinha. 

0 comportamento desta pluma e analisado tomando-se como 
base os dados de temperature da superfTcie da agua, a temperatura do ar 
e a situagao sinotica. A pluma apresentou uma extensao de aproximadamen^ 
te 50 km, com variagoes em seu contorno. Analisando-se a fotografia, ve 
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rifica-se que existe uma wudanga na dlregao da p1un)a> em torno de 13 
km de distancla da fonte, No dia 15/4/75 houve a passagem de uma fren^ 
te fria na regiao; ventos de SW e S pens tatl; ram ate o priroeiro pert£ 
do do dia 17/4 (Tahela 3J, quando ventos de E, NE e N comegaram a s£ 
prar» No instante da passagem do satellte GMT), a diregao da 
pluma indicava ventos de NNE, enquanto para dlstancias malores que 13 
km, a diregao da pluma esta modificada, Indicando que a diregao do yen 
.to tendera para N. 0 efeito de oscllagao da diregao do vento na dire 
gao da pluma fol significativoj desta forma, os valores cal cul ados pe 
las equagoes nao podem ser estendidos para dlstancias malores que 13 
km, devido Ss restrigoes Impostas pela condigao estaclonaria da dire 
gao do vento, razao pela qual restringem-se os comentarios para dista£ 
das menores que 10 km. 

A tempera tura da agua superficial (Tw) na regiao de Ca 
bo Frio no instante da passagem do satellte era de 22,50C, enquanto a 
temperatura do ar (Ta) era de 25,0oc, For Isso a pluma apresenta-se 
multo bem definida e com caracterTsticas de establlidade neutra. 

Caso n9 2 


Com uma extensao em torno de 130 km, a pluma obtida no 
dia 23/6/76 apresentou uma configuragao bastante e^cavel, com contorno 
bastante definido (Figura 5|, Os valores de medidos na fotografla 
foram menores que os previstos pelos metodos propostos (figura 6),sendo 
qiie a categoria. de establlidade prevista pelo metodo de P-G foi neutra 
(D). Os valores de Oy medidos apresentaram caracterTsticas quase neutras 
para dlstancias proximas da fonte e assumiram caracterTsticas estaveis, 
ficando abalxo dos valores previstos pela categorla F de P-G para as 
dlstancias maiores que 5 km. Todos os metodos superestlmaram os valores 
de Oy medidos, em todo o interval o (.Tabela 4). 

Dados de temperature da agua superficial para o perTodo 
de 16/6 a 23/6 (label a 3) forneceram uma temperatura media de 22, 60 C , 
enquanto a temperatura media do ar foi de 22,79 C. No instante da pa£ 
sagem do satelite (11:38 GMT) registrou-se = 22, Oo C e T^ = 23, Qo C, 
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Figura 5 - Fotografia obtida da imagem do MSS 
do LANDSAT ; Caso n9 2. 
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Figura 6 - Valores de dy medidos atraves de fotografia 
em fungao da'^distancia da fonte, cotnparados 
com valores previstos pel os metodos de Pa^ 
quill -Gifford, JUlich e Hogstrbm. 
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alem de se verificar que a elevagao sqlar fgt de 22° e a diregao e In 
tensldade do vento foram de NNE e 8,3 respectlvamente. Dados de 
radiossondageni do Galeao inostraran) a existencia de uma Inversao con) 
base na supar^fTcle qua se estende ate 25Q m de altura, .aclma da qual 
existia uma camada neutra« A mesma situagao fora notada no dia ante 
rior, A diregao do vento sofreu varlagoes no perTodo de 17/6 a 2Q/6,de 
vido a entrada de uma frente que acabou se dlsslpando no dia 19/6. A 
parti r do dia 20/6 houve a predoininancia do anticiclone subtropical 
que se estendeu sobre parte do contlnente* predomlnando em toda reglao 
costelra. A partir do dia 21/6 a diregao do vento passou a ser de NNE, 
com varlagoes noturnas para N, com altas velocldades (Tabela 3). Por 
Isso a pluma apresentou*^se multo bem definida e com uma longa extensao, 

Caso nQ 3 


Neste caso, a fotografla obtida da Iroagemdodia 22/5/77 
(Figura 7) mostra uma pluma de curta extensao (- 20 km) com caracterT^ 
ticas nao muito bem definidas. 0 contorno estabelecldo para ’^sta p1]j 
ma nao representa a realidade, em termos de difusao para os metodos 
existentes’, que cons Ide ram a diregao do vento constante. Este caso es 
ta sendo apresentado para evidenclar os erros que poderao ocorrer nas 
medidas de Oy, quando as condigoes de contorno da pluma nao forem muj^ 
to bem definidas. 

Houve no penodo de 16/5 a 22/S muitas varlagoes na dj_ 
regao do vento (Tabela 3), o que pode ser observado na fotografla apre 
sentada na h‘gura 7, onda a pluma aparece mals bem definida ate a d1^ 
tancia de 5 km da fonte, o que caracteriza urn perTodo de aproximad^ 
mente 1 hora com ventos de SW. Para este mesmo perTodo a temperatura 
da agua superficial apresentou-se maior que a temperatura do ar; no 
Instante da passagem do satellte registrou-se « 23, Qo C e Tg^=*21,0o 
C. No dia 17/5 houve a passagem de uma frente fria na regiao, que oc^ 
sionou varlagoes na diregao do vento. 

I 

Considera-se que os resultados apresentados na Figura 8 
nao podem ser comparados com os prevlstos pelo metodos propostos, ra 
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Figura 7 - Fotografia obtida da imagem do MSS 
do LANDSAT; Caso n9 3. 
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Figura 8 


Valores de ay medidos atraves da fotografia 
em fungao da di stand a da fonte. comparados 
com valores previstos pel os metodos de Pa^ 
quIll-Gifford, OUlich e Hogstrdm. 
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zao pel a qual se deixa de apresentar o calculo das porcentagens de erro 
(label a 4) devtdo ao fato de nao serem representatives de uma situagao 
normal. No entanto, documentam uma situagao que pode proper novos crite 
rios para o calculo de concentragao de poluentes para longas d1 stand as 
da fonte, levando-se em conta as variagoes que ocorrem na dIregSo do 
vento e a sua velocidade, fato que ocorre constantemente no me1o ambien 
te. 


7. Conclusoes 


Anal1sando-se os resultados alcangados neste trabalho, 
tiram-se as seguintes conclusoes: 

A entrada da f rente fria na regiad de Cabo Frio provoca vari£ 
goes na diregao da pi uma, destruindo a sua configuragao inicial 
e formando uma pi uma que,ao longo de sua extensao, apresenta dj_ 
regoes alternantes (Case n9 3), 

“ A difusao sob re o oceano difere apreciavelmente da que ocorre 
sobre o continente, e o comportamento da pi uma e altamente de 
pendente da estab.11 Idade do ar sobre o ocearo, principalmente 
nos baixos nTveis, onde a difusao e governada pel a diferenga de 
temperatura existente entre a agua e o ar. 

- A pi uma mais extensa foi verificada para uma condigao estavel 
(Caso n9 2), alcangando 130 km de comprimento sobre o oceano. 
Para este caso verificou-se que a velocidade do vento foi malor 
que 8 ms""^ e que a diregao do vento permaneceu aproximadamente 
constante em torno de NNE. 

- Embora neste trabalho tenham sido usados os metodos de Pasquill 
-Gifford, JUlicK e Hdgstrdm para comparer os valores de Oy me 
didos atraves das fotografias, cabe ressaltar que esses meto 
dos sao bastante limitados, visto que foram formal izados para 
serem usados sobre o continente e para emissoes de curta dur£ 
gao. Com relagao ao metodo de Pasquill -Gifford, os valores me 
didos nas fotografias assumiram valores que estivaram entre 
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1/3 e 1/2 dos previstos, dependendo da classe de establlldade 
e da d1 stand a da font6, 

0 roetodo de Hdgstrdm demonstra ser eficlente na avallagao de 
Oy para o case quase neutro (Figura 4i, ao contrario do caso 
estavel (Figura 6) que apresentou alta porcentagem de erro (Ta 

bela 31, 

/ 

- 0 metodo de JUlich apresentou porcentagem de erro que varlou eni 
tre 10Q e 300?i, demons trando que sua utilizagao e restrita a 
locals que apresentem rugosidade superficial a1ta*Portanto,quan 
do em alguma situagao estes metodos forem utilizados sobre siu 
perfTcie da.agua, sugere-se que os valores de Oy obtidos fjejam 
devidamente corrigidos. 

Este trabalho demonstrou que a utilizagao do Sensoriament 
to Remoto, atraves de imagens de satelites, para o estudo da difusao a^ 
mosferica e de suma importancia e aplicabil idade pratica, principalmen 
te sobre o oceano, onde as plumas alcangam longas extensoes e pelo fato 
de nao existirem estudos de difusao em larga escala sobre a agua, comp^ 
ravels com aqueles realizados sobre o continente. Acrescente-se o fato 
da excel ente qual idade das fotografias obtidas, devido ao contraste e^ 
tre a pluma e a superfTcie de agua, possibil itando determinar os contor 
no’s da pluma de uma forma contTnua, Esta tecnicasefa de grande utilid^ 
de para estudos futures, permitindo que se estabelegam model os para dj_ 
fe rentes regioes. 
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